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Abstract 

We discuss how it is possible to obtain a reliable approximation 
for the density of states for a system of particles in an anisotropic 
harmonic oscillator potential. A direct application of the result to 
study Bose-Einstein condensation of atomic gases in a potential trap 
can be given. In contrast to a previous study, our method involves 
only analytic calculations. 



Recent experiments on cold alkali gases ^ have given the best evi- 
dence so far for Bose- Einstein condensation (BEC). The experiments demon- 
strate that below a characteristic temperature on the order of a few degrees 
microKelvin several thousand atoms are found in the ground state. Because 
the gases are dilute, as a first approximation interatomic interactions should 
be negligible. However the free boson gas model as described in the stan- 
dard textbooks [Q is not a good theoretical model for these atomic systems 
due to the effects of the confining potential traps used in the experiments. 
To a good approximation the magnetic traps can be modelled by harmonic 
oscillator potentials. There have been several studies of BEC for a system of 
particles in a harmonic oscillator potential ||, |l^, [11 . 



If we consider a system of particles in an external anisotropic harmonic 
oscillator potential characterized by angular frequencies uJi,U2 and u^, then 
the energy levels are (with h = 1) 

111 

^nmana = K + + K + ^)i^2 + K + -)c^3 , (l) 

where /ii,2,3 = 0, 1, 2, . . .. If we describe the system using the grand canoni- 
cal ensemble (see Ref. U2[ for a treatment using the microcanonical ensem- 
ble) , then the thermodynamics may be studied from the grand potential (see 
Pathria in Ref. ||]) 

q=- E In [l - e"^^^"i"2"3-/^)l . (2) 

ni,n2,n3 

The number of particles is 



".1,^2, "3 



It is not possible to evaluate the sums in (^ or (|]) directly in closed form. 
There have been two main approaches to try and obtain analytic approxi- 
mations for such sums. 

One method is to argue that for situations of experimental interest since 
the temperature T is of the order a few microkelvin and Ui is of the order a 
few hundred Hz, we have jSuJi <C 1. It is therefore possible to approximate 
the sums over discrete energy levels in expressions such as (H) and (|]) by 
integrals. A crucial feature in obtaining a reliable approximation is to use an 
appropriate density of states. A beautiful and careful account of the correct 
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way to analyze this problem was given by Grossmann and Holthaus ||^ , and 
used in Ref. 0. These authors noted that the density of states included 
an additional term beyond that found in the unconfined Bose gas. They 
postulated that the density of states was given approximately by 

where f2 = {uJiUJ2UJ3y^^. If only the first term on the right hand side is 
included as in Refs. |^ |[ it leads to an unphysical discontinuous behaviour 
in the specific heat. The second term on the right hand side of @ involves 
a constant 7 which for the isotropic harmonic oscillator is easily seen to be 
3/2, but which Grossmann and Holthaus had to evaluate numerically in the 
anisotropic case. 

A different method for studying and was used by the present 



authors fl^, |TI| . Instead of converting the sums directly into integrals using 
a density of states factor, a contour integral representation was found for 
and d^) with the sums being performed exactly. By deforming the integration 
contour in an appropriate way asymptotic expansions were found for q,N and 
other thermodynamic quantities. These calculations, in contrast with those 
of Refs. [|7|, ^ , were entirely analytic. Given that the approaches of Refs. |]^, |^ 
and [p!0| , pj]] both give good agreement with a direct numerical evaluation of 
the sums (0) and (|), this indicates that by comparing the two sets of results 
for say the particle number, it is possible to extract an analytic result for 7 
in the anisotropic case. As this is rather a roundabout way of proceeding, 
in the present paper we wish to describe a more direct way of evaluating 
the density of states, the result of which leads to an exact equality between 
the approach of Grossmann and Holthaus ||^, |^ and that of our earlier work 

0,0- 

The basic idea involves the use of the asymptotic high-temperature ex- 
pansion of the partition function Z{(3) for the simple harmonic oscillator. 
Generally, if we define Ej\f to be the eigenvalues, labelled by some index N 
{Ejy < Ejy+i), of some self-adjoint differential operator, we may form the 



sum 



POO 

m = E e-"''- = / e-^''p4E)dE , (5) 

with < 6 < £"0 and pd{E) = J^N^iE — E]\r) is the density of states. 
Approximating Pd{E) by a smooth function Ps{E), a refinement of Karamatas 
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theorem states [IT^ that if 



ZiP) = r e-f''' ps{E)dE ~ QP-"^^ + , (6) 

"'^ i=i 

for < (3 < C for some real C > and r^+i < < • • ■ < ri, then for E > b 

k 

Ps{E) = Y.T^,E^'-' + 0{E^^-^~HnE), (7) 

where l/r(rj) is defined to be zero if is a non-positive integer. By virtue 
of this theorem, an asymptotic expansion for the Laplace transform of the 
density of states is carried over into an expansion for the density of states 
itself. Whereas the Laplace transform method can only be used if rj > 0, 
the above theorem extends to all values of r^. An excellent and accessible 
discussion of this theorem with extensive references and applications is the 



book by Baltes and Hilf WA 



The coefficients Cj in (Jg) are known for a wide class of operators on general 



Riemannian manifolds with a variety of boundary conditions [|T5|. However 
in the case of the anisotropic oscillator a direct evaluation of the q is very 
simple. For the eigenvalues E]s[ in we take the energy levels given in 
However, because we are interested in BEC where the ground state plays a 
key role, we single out the ground state for special treatment. If the ground 
state is separated off from the sums in @) and @, then the partition sum 
for determining the density of states is 

oo 3 _-|^ 

Z{13) = e-^("i'^i+"2^2+"^'^^) = n (l - e"'^''') ■ (8) 

ni,n2,?i3=0 i=l 

The case of the harmonic oscillator potential is particularly simple since the 
sum in (|^) involves only geometric series. It is now straightforward to expand 
(H) in powers of t using 



(1 



(The coefficients appearing in the expansion are just the Bernoulli numbers.) 
It is easily found that 

Z{(3) ~ ci/5-=^ + + c,(3-^ + . . . , (10) 
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where 



Ci = {LJ1LJ2UJ3) ^ (11) 

. (12) 



2 \uJ1UJ2 UJ2UJ3 UJ1UJ3, 

1 f ^3 , ^2 , ^1 3 3 3 \ 
^3 = TTT + + + — + — + — . (13 

12 \UJ1UJ2 t^lt^S t^2^3 ^1 ^2 ^3/ 

Using d^) the density of states is found to be 

p{E) ~ ^ciE^ + C2E + C3 . (14) 

Comparison of our result with the ansatz of Grossmann and Holthaus in 
d^) shows that the coefficient of E"^ is correct, and in addition fixes their con- 
stant 7 in terms of the oscillator frequencies. There is an additional constant 
term in the density of states which will make a small correction to the results 
of Refs. [0, 1^ . However, slightly below the critical temperature Tc this con- 
version error is expected to be of the order of or less than the discretization 
error and it may be consistently neglected excluding a temperature interval 
around of the order N~^^^ from the analysis In the special case of 
the isotropic harmonic oscillator with cui = 0^2 = cus = it can be seen that 
our result reduces to ci = Q^^,C2 = and C3 = Q^^. This is also in 

agreement with the results of Ref . . 

Use of the modified form of the density of states (|T^) to discuss BEC 
is basically identical to that found in Refs. 0, |^, and there is no need to 



repeat the analysis here. The coefficients ([TlH lq) also enter the asymptotic 



expansions of Refs. |]T0|, |ri| , and it is now straightforward to establish a direct 



connection between the results of Refs. § and those of Refs. [p^ |TT 



After this work was completed we became aware of another systematic 



treatment of Bose-Einstein condensation in harmonic oscillator traps ||T6 
This reference uses the Euler-Maclaurin summation technique to analyze 
the harmonic oscillator sums. Detailed analysis is given for the isotropic 



harmonic oscillator which appears to agree with our earlier work |T^, |TT| . A 
result for the particle number is given for the case of the anisotropic oscillator 
which is also in agreement with the results of our evaluation |TI|] , as well as 



the result found by using the first two terms in the density of states approach 
given in the present paper. 

In conclusion, we have shown how to calculate the density of states factor 
for a system of particles in an anisotropic harmonic oscillator potential. This 
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means that an analysis similar to that of Grossmann and Holthaus [0, |^ can 
be performed in an entirely analytic manner without any resort to numer- 
ical computation. The method we have outlined can be extended to many 
other situations of interest, and we will present a more complete discussion 
elsewhere. 

Acknowledgements 

D. J. T. would like to thank J. D. Smith for helpful discussions. 

References 

[|] kirstenOtphlOO . physik . uni-leipzig . de. 

[I] d . j . tomsSnewcastle . ac . uk. 

[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and 
E. A. Cornell, Science 269, 198 (1995). 

[2] C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. Hulet, Phys. Rev. 
Lett. 75, 1687 (1995). 

[3] K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Dur- 
fee, D. M. Kurn, and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995). 

[4] See for example F. London, Superfluids II, (John Wiley, New York, 
1954), K. Huang, Statistical Mechanics (John Wiley, New York, 1987) 
or R. K. Pathria, Statistical Mechanics (Pergammon, Oxford, 1972). 

[5] S. R. de Groot, G. J. Hooyman, and C. A. ten Seldam, Proc. Roy. Soc. 
(London) A 203, 266 (1950). 

[6] V. Bagnato, D. E. Pritchard, and D. Kleppner, Phys. Rev. A 35, 4354 
(1987). 

[7] S. Grossmann and M. Holthaus, Phys. Lett. A 208, 188 (1995). 

[8] S. Grossmann and M. Holthaus, Zeit. fiir Naturforschung 50 a, 921 
(1995). 

[9] H. Haugerud and F. Ravndal (unpublished), |cond-mat / 950904T . 



6 



[10] K. Kirsten and D. J. Toms, Bose-Einstein condensation of atomic gases 
in a harmonic oscillator confining potential trap, preprint, Feb. 1996, 
|cond-mat/ 96040311 - 

[11] K. Kirsten and D. J. Toms, Bose-Einstein condensation of atomic gases 
in general harmonic oscillator confining potential traps, preprint, March 
1996, to appear in Phys. Rev. A. 

[12] S. Grossmann and M. Holthaus, Microcanonical fluctuations of a Bose 
system's ground state occupation number, preprint, Feb. 1996. 

[13] F. H. Brownell, Paciflc J. Math. 5, 483 (1955); J. Math. Mech. 6, 119 
(1957). 

[14] H. P. Baltes and E. R. Hilf, Spectra of Finite Systems (Bibliographisches 
Institut Mannheim, Zurich, 1976). 

[15] S. Minakshisundaram, Can. J. Math. 1, 320 (1949); S. Minakshisun- 
daram and A. Pleijel, Can. J. Math. 1, 242 (1949); B. S. DeWitt, Dy- 
namical Theory of Groups and Fields (Gordon and Breach, New York, 
1965); P. B. Gilkey, Invariance theory, the heat equation and the Atiyah- 
Singer index theorem, 2nd. Edn. (CTC Press, Boca Raton, 1995). 

[16] H. Haugerud, T. Haugset, and F. Ravndal, A more accurate analysis 
of Bose-Einstein condensation in harmonic traps, preprint. May 1996, 
|cond-mat/ 9605 log . 



7 



